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Abstract 

We study the squark loop correction to W^H^ associated hadroproduction via 
gluon-gluon fusion within the minimal supersymmetric extension of the standard 
model. We list full analytic results and quantitatively analyze the resulting shift 
in the cross section at the CERN Large Hadron Collider assuming a supergravity- 
inspired scenario. 
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1 Introduction 



The search for Higgs bosons will be among the prime tasks of the CERN Large Hadron 
Collider (LHC) JT|. While the standard model (SM) contains one complex Higgs doublet, 
from which one neutral CP-even Higgs boson emerges in the physical particle spectrum 
after the electroweak symmetry breaking, the Higgs sector of the minimal supersymmetric 
extension of the SM (MSSM) consists of a two-Higgs-doublet model (2HDM) and accom- 
modates five physical Higgs bosons: the neutral CP-even h° and H° bosons, the neutral 
CP-odd A boson, and the charged if ± -boson pair. At the tree level, the MSSM Higgs 
sector has two free parameters, which are usually taken to be the mass rriA of the A 
boson and the ratio tan/5 = Vz/vx of the vacuum expectation values of the two Higgs 
doublets. 

The discovery of the if ± bosons would rule out the SM and, at the same time, give 
strong support to the MSSM. The main strategies for the if ± -boson search at the LHC 
were summarized in Refs. Depending on the if ± -boson mass m#, the dominant 

mechanism of single if ± -boson hadroproduction are gg, qq — > ti followed by t — > bH + U, 
gb — > tH + ||, gg — > tbH + ||J], and qb — > q'bH + || together with their charge-conjugate 
counterparts. The hadroproduction of H + H~ pairs proceeds at tree the level via qq 
annihilation, qq — > H + H~ , where q = u,d,s,c |J, and b 0, and at the one-loop level via 
gg fusion, gg — > H + H~, which is mediated by quark ]7],|3||| and squark loops J7|,^|. The 
suppression of the gg-fusion cross section by two powers of the strong-coupling constant a s 
relative to the one of qq annihilation is partly compensated at multi-TeV hadron colliders 
by the overwhelming gluon luminosity. 

An interesting alternative is to produce H bosons in association with W T bosons, 
so that the leptonic decays of the latter may serve as a trigger for the if ± -boson search. 
The dominant partonic subprocesses of W ± H T associated production are bb — > W^H^ 
at the tree level and gg — *> W ± H T at one loop, which were investigated for vanishing 
bottom-quark mass m\, and small values of tan/3 (0.3 < tan/? < 2.3) in Ref. 0] and 
recently, without these restrictions, in Refs. PJTT|]. A careful signal- versus-background 
analysis, based on the analytic results of Ref. 0, was recently reported in Ref. [0. So 



far, only the quark loop contribution to gg — > H /± if =F was considered [Q,|T^,|n|]. The 
purpose of this paper is to provide, in analytic form, the supersymmetric contribution 
to this partonic subprocess, which is induced by virtual squarks through the Feynman 
diagrams depicted in Fig. |1|. Furthermore, we wish to quantitatively study its influence 
on the cross section of the inclusive reaction pp —> W ± H T + X at the LHC. We recall 
that, in the case of pp — > H + H~~ + X, the supersymmetric correction to the gg-fusion 
cross section can be as large as +50% 0. A priori, one expects to encounter a similar 
situation for pp — > W ± H T + X. 

In order to reduce the number of unknown supersymmetric input parameters, we 
adopt a scenario where the MSSM is embedded in a grand unified theory (GUT) involving 
supergravity (SUGRA) |13|]. The MSSM thus constrained is characterized by the following 
parameters at the GUT scale, which come in addition to tan (3 and rriA'- the universal scalar 
mass mo, the universal gaugino mass m^, the trilinear Higgs-sfermion coupling A, the 
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bilinear Higgs coupling B, and the Higgs-higgsino mass parameter \i. Notice that ttla is 
then not an independent parameter anymore, but it is fixed through the renormalization 
group equation. The number of parameters can be further reduced by making additional 
assumptions. Unification of the r-lepton and 6-quark Yukawa couplings at the GUT scale 
leads to a correlation between m t and tan (3. Furthermore, if the electroweak symmetry 
is broken radiatively, then B and \x are determined up to the sign of \i. Finally, it turns 
out that the MSSM parameters are nearly independent of the value of A, as long as 
\A\ ^500 GeV at the GUT scale. 

This paper is organized as follows. In Sec. §, we list the helicity amplitudes of the 
partonic subprocess gg — > W~H + involving virtual squarks. In Sec. |3|, we present quanti- 
tative predictions for the inclusive cross section of pp — > W ± H T + X at the LHC adopting 
the SUGRA-inspired MSSM. Sec. f| contains our conclusions. 



2 Analytic Results 

In this section, we express the gg — > W~H + helicity amplitudes involving one closed 
squark loop in terms of the standard scalar two-, three-, and four-point functions, 



pl,ml,ml 



5o(; 

Co (pi, (P2-Pi) 2 ,pI 



d D q 1 

in 2 (q 2 — ml + ie) [(q + pi) 2 — m\ + ie] 



d D q 



in 2 (q 2 — ml + ie) [(q + p\) 2 — m\ + ie] [(q + p<i) 2 — m\ + ie] 
Do (pi, {P2~Pif, {P3~P2) 2 ,pI,pI, {p3-pi) 2 ,m 2 ,ml,ml,mfj 
d D q 1 



in 2 (q 2 — ml + ie) [(q + p\) 2 — m\ + ie] [(q + p 2 ) 2 — m 2 , + ie] [(q + p^) 2 — m\ + ie] ' 

(1) 

where D is the space-time dimensionality. The Bq function is ultraviolet (UV) diver- 
gent in the physical limit D — » 4, while the Co and Dq functions are UV finite in 
this limit. We evaluate the Bq, Co, and Do functions numerically with the aid of 



the program package FF [fLi}] . To simplify to notation, we introduce the abbreviations 



C°j k (c) = C (a, b, c, mf, m 2 , m|) and Dfj^(e, f) = D (a, b, c, d, e, /, m;, nr;. mi. in 



We work in the MSSM adopting the Feynman rules from Ref. |jl5| . For each quark 



flavor q there is a corresponding squark flavor q, which comes in two mass eigenstates 
z = l,2. In the following, up- and down-type squark flavors are generically denoted by t 
and b, respectively. The masses m^ of the squarks and their trilinear couplings g w ^.i., 
9h°qiqji dHOqigp and g H +i~ h . to the W~ , h°, H°, and H + bosons are defined in Eqs. (A. 5) 
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and (A. 9) of Ref. [|16| and in Eq. (A. 2) of Ref. |7J, respectively. 1 Furthermore, we have 

cos(a — (3) 



9W~H+h° 
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sm(a - (3) 

9W-H+H0 = - , [ZJ 

where a is the mixing angle that rotates the weak CP-even Higgs eigenstates into the 
mass eigenstates h° and H°. 

Calling the four-momenta of the two gluons and the W boson p a , pb, and pw, respec- 
tively, we define the partonic Mandelstam variables as s = (p a + Pb) 2 , t = (p a — pw) 2 , an d 
u = (pb — pw) 2 - Furthermore, we introduce the following short-hand notations: w = m^, 
h = m 2 H , d = t — u, ti = t — w, t 2 = t — h, Ui = u — w, u 2 = u — h, N = tu — wh, 
A = s 2 +w 2 +h 2 — 2(sw+wh+hs), and q = m~ — m| . We label the helicity states of the two 

gluons and the W boson in the partonic center-of-mass (cm.) frame by A a = —1/2, 1/2, 
X b = -1/2, 1/2, and \ w = -1, 0, 1. 

The relevant Feynman diagrams are depicted in Fig. [I]. I n analogy to the quark case, 
we refer to the diagrams involving a neutral Higgs boson in the s channel as triangle 
diagrams. In contrast to the quark case, the diagrams involving the A boson add up to 
zero. The residual diagrams are regarded to be of the box type. The helicity amplitudes 
of the squark triangle contribution read 

Mi Xb0 = aVx(i + x a \b) E E ( 9w - 2 H+ + 9h °"t + 9w T H ° 9m % ) 

b 5i\ s ~ m h° + im h°^ h o s - mjjo + im H oT H o J 



x 



M$ aXb±1 = 0, (3) 

where T^o and T^o are the total decay widths of the h° and H° bosons, respectively. In 
this case, the W boson can only be longitudinally polarized because it couples to two Higgs 
bosons. As for the squark box contribution, all twelve helicity amplitudes contribute. Due 
to BoseP] and weak-isospin symmetry, they are related by 

Ml XbXw (t,u,ml,ml) = (-l) Xw M D XbXaXw (n,t,m|,m|) , 

M° a x b x w (^^ m l m i) = ~M-x a -x b -x w (t,u,ml m?) , (4) 



respectively. Keeping Xw = ±1 generic, we thus only need to specify four expressions. 
These read: 



m d ++q = —f=, e T,9 W -u bj 9 H+ ub 3 [K+ + (t 

SVA (Lb) M 



U 



1 In Ref. Q, mfo and g w -i.i. are called Mg a and Vff aDb /g, respectively. 

2 Notice that the interchange of t and u also affects the representation of the VF-boson polarization 
four- vector through its dependence on the angle between the three-momenta of gluon a and the W boson. 
This explains the sign factor in the first line of Eq. (^) , which is not expected from pure Bose symmetry. 
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where denotes the sum over squark generations and 
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Notice that the UV divergences of F}_ and F 2 _ cancel in the expression for A4 
Eq. (|). 

The differential cross section of the partonic subprocess gg — * W~H + is then given by 
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where a s (/x r ) is the strong-coupling constant at renormalization scale /i r , Gp is Fermi's 
constant, and \ and -M\ a \ \ are the helicity amplitudes of the quark triangle and 

box contributions, which may be found in Eqs. (1) and (3) of Ref. [llj, respectively. The 
relative minus signs between the quark and squark terms in Eq. (|7|) compensate for the fact 
that the Feynman rules underlying Ref. [11] differ from those adopted here. Due to Bose 
symmetry, the cross section da/dt of gg — > W~H + is symmetric in t and u. Due to charge- 
conjugation invariance, it coincides with the one of gg — > W + H~ , so that the cross section 
da/dt of gg — > W ± H T emerges from the right-hand side of Eq. (0) by multiplication with 
two. The kinematics of the inclusive reaction AB — > W ± H T + X, where A and B are 
colliding hadrons, is described in Sec. II of Ref. |@. Its double-differential cross section 
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d 2 a / dy dpx , where y and pr are the rapidity and transverse momentum of the W boson 
in the cm. system of the hadronic collision, may be evaluated from Eq. (2.1) of Ref. 0. 



3 Phenomenological Implications 

We are now in a position to explore the phenomenological implications of our results. 
The SM input parameters for our numerical analysis are taken to be Gp = 1.16639 • 
10~ 5 GeV~ 2 , m w = 80.419 GeV, m z = 91.1882 GeV, m t = 174.3 GeV , and m b = 4.6 GeV 
|17|| . We adopt the lowest-order set CTEQ5L of parton density functions for the 



proton. We evaluate a s (/x r ) from the lowest-order formula |l7j with rif = 5 quark flavors 
and asymptotic scale parameter Aqq D = 146 MeV [TH|. We identify the renormalization 
and factorization scales with the W ± H T invariant mass s. For our purposes, it is useful 
to replace by m#, the mass of the H ± bosons to be produced, in the set of MSSM 
input parameters. We vary tan/? and m# in the ranges 2.5 < tan/3 < 38 ~ m t /m b 
and 120 GeV < m H < 600 GeV, respectively. As for the GUT parameters, we choose 
rrii/2 = 150 GeV, A = 0, and n < 0, and tune m so as to be consistent with the 
desired value of ran- All other MSSM parameters are then determined according to the 
SUGRA-inspired scenario as implemented in the program package SUSPECT [[Tj|. We 
do not impose the unification of the r-lepton and 6-quark Yukawa couplings at the GUT 
scale, which would just constrain the allowed tan/5 range without any visible effect on 
the results for these values of tan/3. We exclude solutions which do not comply with 
the present experimental lower mass bounds of the sfermions, charginos, neutralinos, and 
Higgs bosons [pOR . 

We now study pp — > W ± H T + X at the LHC, with cm. energy y/S = U TeV. The 
fully integrated cross section is considered as a function of ran for tan/3 = 3, 10, and 30 
in Fig. |^(a) and as a function of tan/3 for van = 150, 300, and 600 GeV in Fig. ||(b). 
The combined gg-fusion contribution due to quarks and squarks (solid lines) is compared 
with the one due to quarks only (dotted lines) @ . For reference, also the ^-annihilation 
contribution (dashed lines) is shown gj. We note that the SUGRA-inspired MSSM with 
our choice of input parameters does not permit tan (3 and mjj to be simultaneously small, 
due to the experimental selectron mass lower bound PD |. This explains why the curves 



for tan/3 = 3 in Fig. 0(a) only start at ran « 260 GeV and those for ran = 150 GeV 
in Fig. |2|(b) at tan/3 w 9. For large m#, the experimental lower bound [^{J enforces 
tan /3 :> 2.5. On the other hand, the experimental lower bounds on the chargino and 
neutralino masses induce an upper limit on tan (3, which depends on ran- We observe 
from Figs. 0(a) and (b) that the supersymmetric correction to the gg-fusion cross section 
can be of either sign and have a magnitude of order 10%. It exceeds +10% for small tan f3 
and large m#, while it almost reaches —10% for medium tan/3 and small or medium mjj. 
On the other hand, it is generally small for large tan/3. We recall that, in the case of 
pp — > H + H~ + X, the supersymmetric correction to the gg-fusion cross section can be 
as large as +50% 0. As explained in Ref. @, the dip in the ttlh dependence of the 
gg-fusion cross section located about ran = mt [see Fig. H(a)] arises from resonating top- 
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quark propagators in the quark box diagrams. Furthermore, the minima of the curves 
in Fig. H](b) close to tan/3 ~ ^mt/rrih ~ 6 may be understood by observing that the 

average strength of the H~bt coupling, which is proportional to J mf cot 2 (3 + m 2 tan 2 /3, 
is then minimal JF]. As is the quark case @, the squark triangle and box contributions 
are similar in size and destructively interfere with each other, so that their superposition 
is much smaller than each of them separately. As in the case of pp — > H + H~ + X J7], 
the bulk of the squark contribution comes from the stop and sbottom squarks, while the 
contributions from the first- and second-generation squarks is greatly suppressed because 
their couplings to the Higgs bosons are significantly smaller than those of the third- 
generation squarks and their masses are generally larger than those of lightest stop and 
sbottom squarks, t\ and b\. We conclude that the suppression of the gg-fusion cross 
section relative to the one of bb annihilation remains after the inclusion of the squark loop 
contributions. 

It is interesting to find out how the kinematic behavior of the gg-fusion cross section 
is affected by the supersymmetric correction. To that end, we study in Figs. |3](a) and (b) 
the distributions in the ly-boson transverse momentum p? and rapidity y, respectively, 
for tan/3 = 3, 10, 30, and m H = 300 GeV. While the y distribution does not exhibit any 
striking features, we observe that the supersymmetric correction leads to an increase of 
the pt distribution by more than 50% at large pt for medium to large tan f3. This can 
be traced to the presence of absorptive parts in the squark loop contribution. In fact, if 

Pt > J A (4m?, m 2 ^, m 2 ^ / (ArriqJ, then s > 2mg., so that pairs of real qi squarks can be 

produced. 

For a comparison with future experimental data, the ^-annihilation and gg-fusion 
channels should be combined. From Fig. 0(a), we read off that the total cross section 
of pp — > W ± H Zf + X at the LHC is predicted to be approximately 500 fb (20 fb) in the 
considered MSSM scenario if tan (3 = 30 and m H = 150 GeV (tan (3 = 3 and m H = 
300 GeV). If we assume the integrated luminosity per year to be at its design value of 
L = 100 fb -1 for each of the two LHC experiments, ATLAS and CMS, then this translates 
into about 100.000 (4.000) signal events per year. 

4 Conclusions 

We calculated the squark loop contribution to the partonic subprocess gg — > W ± H T 
within the MSSM, and analyzed its impact on the inclusive cross section of pp — > W ± H^ + 
X and its distributions in transverse momentum and rapidity at the LHC adopting a 
SUGRA-inspired scenario. Its inclusion may increase or decrease the integrated gg-fusion 
cross section by up to 10%, depending on the values tan/3 and m#. However, bb anni- 
hilation remains the dominant mechanism of W^H^ associated hadropro duct ion at the 
LHC. Should the MSSM be realized in nature, then the W ± H T channel will provide a 
copious source of charged Higgs bosons at the LHC, with an annual yield of up to 100.000 
signal events. 
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Figure 1: One-loop Feynman diagrams for gg — > W H + due to virtual squarks in the 
MSSM. 
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Figure 2: Total cross sections a (in fb) of pp — > W ± H T + X via bb annihilation (dashed 
lines) and gg fusion (solid lines) at the LHC (a) as functions of ran for tan/3 = 3, 10, and 
30; and (b) as functions of tan/3 for m# = 150, 300, and 600 GeV. For comparison, also 
the quark loop contribution to gg fusion (dotted lines) is shown. 
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Figure 3: (a) pr distributions da/dpr (in fb/GeV) and (b) y distributions da/dy (in fb) 
of pp — > W ± H Zf + X via bb annihilation (dashed lines) and gg fusion (solid lines) at the 
LHC for tan/3 = 3, 10, 30, and m# = 300 GeV. For comparison, also the quark loop 
contribution to gg fusion (dotted lines) is shown. 
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